Experimentally friendly bounds on non-Gaussian entanglement from second moments 
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We demonstrate that the entanglement in a class of two-mode non-Gaussian states obtained by subtracting 
photons from Gaussian twin beams can be bounded from above and from below by functionals of the second 
moments only. Knowledge of the covariance matrix thus suffices for an entanglement quantification with ap- 
preciable precision. The absolute error in the entanglement estimation scales with the non-Gaussianity of the 
considered states. 
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I. INTRODUCTION 

Interfaces between light and matter are key building blocks 
of a future quantum web, a global secure communication net- 
work where the manipulation and transmission of information 
are regulated by the quantum laws 11J]. The transfer of quan- 
tum states and the distribution of correlations across the in- 
terfaces are enabled by the common mathematical language 
of canonically conjugate observables with continuous spec- 
trum, such as the quadratures of light and the collective spin 
components of atomic ensembles. It is thus very fascinating 
to witness how the second generation of quantum informa- 
tion research is focusing more and more, from both theoretical 
and experimental viewpoints, on the characterization of con- 
tinuous variable (CV) entanglement and its applications for 
communication, computation and metrology ||2|] • 

Out of the infinite-dimensional Hilbert space of CV states, 
a special class of states has played a prominent role in recent 
years: Gaussian states. Their mathematical treatment is ad- 
vantageous thanks to a compact formalism based on symplec- 
tic analysis, and a very accurate degree of control is reached 
in their experimental realizations with light and matter JH. 
There are however many tasks which are impossible by using 
only the Gaussian states and operations toolbox (e.g., entan- 
glement distillation 14J] and universal quantum computation 
15D), and many other tasks which can be sharply improved 
by suitably resorting to some non-Gaussianity (e.g. CV tele- 
portation (a] and loss estimation |7[]). These premises have 
spurred astonishing progresses in the experimental engineer- 
ing of non-Gaussian states, such as Fock states and states 
obtained by deGaussifying Gaussian resources via addition 
and/or subtraction of single photons [8]. It has been in par- 
ticular verified experimentally that a photon subtraction from 
a two-mode squeezed entangled Gaussian state leads to an en- 
hancement of the entanglement (at fixed squeezing) and 
it is known that such resource might be used for a more per- 
formant quantum teleportation of classical and nonclassical 
states (a primitive of the quantum internet) ll^ [Tol[Tlll and for 
yet-to-be-achieved demonstrations of loophole-free Bell tests 
of nonlocality flZl . The bottleneck for unleashing the power 
of non-Gaussian CV quantum technology has remained the 
quantitative characterization of entanglement in states which 
deviate from Gaussianity, a crucial step to evaluate and con- 
trol their usefulness for applications and the bona-fideness of 



their preparation. While for Gaussian states all the informa- 
tion is encoded in the second moments of the canonical oper- 
ators (collected in the covariance matrix), for any other state 
an infinite hierarchy of moments is in principle needed for an 
exact entanglement quantification 111 311 . This translates, in ex- 
perimental terms, into the demand for a complete state tomog- 
raphy [14], a process which is time- and resource-consuming 
especially for two or more modes. 

In this paper we provide an advance in the characterization 
of non-Gaussian entanglement which reduces the complexity 
of its experimental determination exactly to the same level of 
Gaussian states. This cannot be possible for any generic non- 
Gaussian state: hence here we focus on the important class of 
pho ton-subtracted states (PSS) under quite realistic conditions 
lllOillin . which represent the preferred resources for most cur- 
rent and future applications requiring non-Gaussianity. Com- 
bining results from the extremality of Gaussian states lfl5ll 
with an analysis of quadrature correlations, we derive ana- 
lytical lower and upper bounds which individuate the entropy 
of entanglement of pure two-mode PSS (obtained from Gaus- 
sian twin beams by conditional subtraction of k photons per 
beam) within a very narrow band, with relative error vanish- 
ing with increasing entanglement. The bandwidth is linked 
to the degree of non-Gaussianity lfl6ll of the analyzed PSS. 
The results are extended to bound the entanglement of forma- 
tion of a class of mixed PSS obtained by means of realistic 
"on/off" type photodetectors. Crucially, all these bounds are 
only functions of the second moments of the non-Gaussian 
states. A novel method for the fast and reliable reconstruc- 
tion of the complete covariance matrix of optical two-mode 
CV states has been recently demonstrated lfl7l [l8tl : our result 
proves that it can be applied to real deGaussified PSS as well 
and it is enough to quantitatively estimate entanglement with 
surprisingly high accuracy. 



II. PRELIMINARIES 

We deal with a CV system of N = 2 bosonic modes, asso- 
ciated to an infinite-dimensional Hilbert space Ti = F% <g> F^ 
[31]. Here Fi is the Fock space of each individual mode i, de- 
scribed by the ladder operators fij , d\ satisfying the canonical 
commutation relations [di, at] = Sij. We can collect the field 
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quadrature operators into the vector X = {qi,Pi,q2,P2}, 
with fa — a,i + a\ and pi = fij — a\. The exact descrip- 
tion of a generic CV (non-Gaussian) state requires arbitrary- 
order moments of the canonical operators. In general, the first 
moments X = ((Xi), (Xi), (X 2 ), (X 2 )) can be adjusted by 
local displacements without affecting entanglement: they will 
be set to zero without loss of generality. Two-mode Gaussian 
states are henceforth completely specified by the 4x4 real 
symmetric covariance matrix (CM) <x of the second moments 
a-ij = (X l X J + XjXi) /2 - (Xi) (Xj) |H. The CM elements 
for an arbitrary two-mode CV state can be efficiently recon- 
structed in the lab by means of homodyne detections lll7i[l^l . 
We recall that any physical two-mode CM 
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can be converted by local unitaries into a standard form [2C 
where on — diag{ai, a^} and 7 = diag{7ir, j p } with a, > 
1,7*>|7pI>0. 

Entanglement of a pure bipartite state \tp) is universally 
quantified by the Von Neumann entropy of the reduced den- 
sity matrix of each subsystem (entropy of entanglement E) 
Pill . The optimal convex-roof extension of E for a mixed 
state g defines the entanglement of formation l2lll Ep(g) = 
infjp. ^.j J^iPiEfyi), where the infimumruns over all pure- 
state decompositions of g = X)i-Pi|V , i)(V'i|- Another popu- 
lar measure of entanglement for pure and mixed states, whose 
computation is in general simpler, is the logarithmic negativity 
JUIll] E N = log [tr\g Ti |] IH, where the partial transposi- 
tion g Ti is obtained from g by transposing the degrees of free- 
dom of one subsystem only. For pure states En > E. An im- 
portant result in CV entanglement theory is the 'extremality' 
of Gaussian states: for any CV state g with second moments 
given by cr, the corresponding Gaussian state g G defined by 
the same CM has smaller entanglement (quantified by a con- 
tinuous and strongly superadditive measure) than g IU5I1 . It is 
also known that for arbitrary two-mode Gaussian states Ep is 
computable, additive, and strongly superadditive I24ll25ll . 
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where 2 Fi denotes the Gauss hypergeometric function. We 
define an effective squeezing parameter r such that 

TX = z = tanh r. 

Thus |V> fc ) =A r fe " 1/2a i a 2l^o(r)), where 



A4 = E 
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n—k 

= fc! 2 2 -Fi (k + l,k+ 1; l;tanh 2 r) cosrT 2 r tanh 2fe r . 

The logarithmic negativity of the PSS states for arbitrary k 
can be computed in simple closed form via the formula 

EN(ipk) — 2 log[5^?Li, ci^] and reads in our notation 



ry k 
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COVARIANCE MATRIX AND ENTANGLEMENT OF 
PURE PHOTON SUBTRACTED STATES 



E N (i> k ) = - log[(l - z) 2fe + 2 2 Fi (k + 1, k + 1; 1; z 2 )] . 

(2) 

Notice that E^(ipk) increases with k: iteration of the photon 
subtraction process further enhances the entanglement com- 
pared to the original Gaussian instance. We aim at an esti- 
mate of the entropy of entanglement of PSS for any k, and 
specifically at accurate bounds on this universal entanglement 
quantifier which can be measured experimentally with high 
efficiency. We will now derive the CM of PSS states in closed 
form and show that it contains enough information, straight- 
forwardly accessible in the lab, for the desired task. 

Denoting by O an observable involving normally ordered 
combinations of ladder operators on modes 1 and 2, we ob- 
serve that (ip k \d\ip k ) = Afk 1 (Ma\ k a 2 k 6d l ld^o), that 
is, we can evaluate expectation values of relevant operators 
in terms of normally ordered higher moments computed on 
the unperturbed Gaussian twin beams. Recalling that, for 
\ijjo(r)), the normally ordered characteristic function reads [0] 

Xq (cti, a*, Q!2, a 2 ) = exp{sinhr[(aiQ;2 + Oi\a 2 ) coshr — 
(|ai| 2 + |c*2| 2 ) sinhr]}, we can write the (normally ordered) 
moment generating formula: 



The starting point for the definition of "ideal" PSS is a 
pure two-mode squeezed Gaussian state (or 'twin beam' in the 
optical language), |-0 o ( r o)) = J2 n ^ n VT ~ A 2 |n, n), where 
A = tanhro and the positive ro is the squeezing degree. 
The CM <T for this state is already in standard form with 
cli = a 2 — cosh2ro, j x — ~7p = sinh2fo. The beam 
1 (2) of \ipo) is let to interfere, via a beam splitter with 
transmittivity T (preferably T close to unity), with a vacuum 
mode 1' (2')- The output is a four-mode Gaussian state of 
modes 1,1', 2, 2'. A photon-number-resolving detection of 
exactly k photons in each of the two beams 1' and 2', con- 
ditionally projects the state of modes 1,2 into a pure sym- 
metric non-Gaussian state 02611 . given in the Fock basis by 



Mt = ^ly(aiT(&tr(a 2 )% 



(3) 
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The formula Eq. © can be readily applied to compute the sec- 
ond moments of the canonical operators on our non-Gaussian 
states. After some algebra, the CM <Jk of PSS of the form 
IV'fc) turns out to be that of a symmetric two-mode squeezed 
thermal state, automatically in standard form, with a\ = 
oW = M k -\M k k > + k 2 , k + M^ k + 2M*+# + M k k] k k ), 
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FIG. 1: (color online). Upper (solid) and lower (dotted) bounds on the entanglement of PSS obtained from second moments, plotted as 
functions of the squeezing r. The entanglement of formation is individuated within the narrow shaded region between the bounds. Panels 
[(a)-(c)] depict ideal PSS |f/>fc)> whose logarithmic negativity is shown as well for reference (thin dashed line). Panel (d) depicts a more realistic 
PSS f?fe =1 modeled as a binomial mixture of the Gaussian twin beam and the first four pure PSS, £>s = i = Ysk=oO- — p) n ~ k p k (^) |y>fc)(?/>fc|, 
with n = p" 1 = 4; notice how only slight modifications occur to the bounds compared to the corresponding ideal case (a). 



^; fc it/c 1 +<^f)- Explicitly: 

x [2(k + l) 2 sinh 2 r 2 Fi(-fc,-fc;2;tanh 2 r) 

+ 2 ^ 1 (-/c,-fc;l;tanh 2 r)]; (4) 
(fe) _ 2 tanhrQ + l) 2 Fi(fc + l,fc + 2;l; tanh 2 r) 
2 Fi(fc + l,A; + l;l;tanh 2 r) 

We will now show how to extract from the CM elements a 
lower and an upper bound on E(ipk)- 

Lower bound. — The extremality of Gaussian states II 1 311 en- 
tails that S 
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where is the mixed Gaussian state with CM er fc . Ex- 
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Upper bound. — The logarithmic negativity of pure PSS is al- 
ready an upper bound for E, however it is defined in terms 
of all the moments of the state, and its experimental determi- 
nation (in real conditions requiring a complete state tomog- 
raphy) becomes rather demanding even when the states take 
very special forms The strength of our investigation is to 
derive a slightly looser upper bound but which is a function 
of the second moments only of the PSS, i.e. of cr k - We sim- 
ply observe that for k = 0, namely for Gaussian twin beams, 
En(iPo) = arcsinh^ ']. The standard-form parameter 



quantifies the maximum quadrature correlation between the 
two modes j28ll . It is tempting to postulate that, for PSS 
with an arbitrary degree k of deGaussification, a function of 
jW ma y yjeid an overestimate of the actual entanglement 
(which, we remark, is nontrivially encoded in higher-order 
correlations too). We can turn this blurry bit of intuition into 
the following 
Theorem 1. For all k > 0, 



E% p = log[l + 2 7 «] > E N (^ k ) > E(<p k ) . 



(7) 



Proof. The rightmost inequality holds by definition. Here we 
sketch the (quite technical) proof of the leftmost one, which 
is one of the main results of this paper. The simple cases 
k = 0, 1 can be proven by inspection, hence we specify here 
to arbitrary k > 2. Using Eq. (|2]l and Eq. (0, and expo- 
nentiating both sides of the inequality, the problem reduces to 
proving that F^(z) = 2 F\ (k + 1, k + 1; 1; z 2 ) + Az(k + 
1) 2-F1 (fc + 2, k + 1; 1; z 2 ) - (1 - z)- 2 *^ 2 > where we 
recall that tanhr = z £ (Q, 1). We can write F^ k '(z) as 
a power series, F^(z) = J2m=i fm 1 z m , where fm = 
[(1 - (-1)-) (2k + m + |) + 1] (i(- 1 +(-i fe ) m )+ fe +^) 2 - 

< Vj > 1, 



( 2 UTi +1 )' We observe that fS-i ^ and ft* 

and moreover > —f%f V? > k + 1. This entails that 

by truncating the power series at m = 2k + 2 we discard a 

positive remainder: F^(z) > F^{z) = Y^=if™ ] z m - 
Let us now define a parametric class of hypergeometric sums, 

Q (k) _ ^2fe+2 (2k+l-m+l\ f {k) 



f(*0 



)/m ■ By means of Zeilberger's 



4 



algorithm JH] one can verify that Sf > Vfc > 2, I > 0. We 
will now show that p( k) (z) > s[ k) z 2k+2 > to conclude 
the proof. We take the ratio (z) = [F (k *> (z)]/[S { k) z 2k+2 } 
and expand it in power series around z = 1~: R( k '(z) = 

1 + E^i(-l)'[Sf V4 fc) ]0 - !)'■ But trivially (z - 1)' = 
(— 1) (1 — zf, and being z < 1 the alternating sign is can- 
celled to yield R (k *> (z) > 1. □ 
We have shown that, quite remarkably, the simple measure- 
ment of the CM of a pure PSS enables to pin down the en- 
tropy of entanglement quantitatively within analytical a priori 
bounds. In fact, one can appreciate how close the lower and 
upper bounds (both functions of the second moments only) 
are to each other for various values of k in Fig. []"[(«)-(£)]. A 
crucial fact is that the absolute error A fe = (E^ p - E l k ow ) /2 on 
the entanglement quantification asymptotically saturates (for 
r — > oo ) to a constant value 

Ar x =5[log(4 + 8fc)-l], 

being Afe < A™ ax for any finite r. Accordingly, since the 
actual value of the entanglement (measured by the average 
between lower and upper bound) diverges linearly with the 
squeezing r, the relative error S k = (E^ p - E k ow )/(E k xp + 
E l k ° w ) on the estimate of E(ipk) from the CM vanishes for 
r ^> 0, rendering our method rigorously accurate. We no- 
tice that, in general, the error Afe increases with k, although it 
stays of the order of few units - on a scale ranging to infinity 
- even for big k (e.g. Afe < 4 for up to k — 1000 photon sub- 
tractions per beam), thus scarcely affecting the quality of the 
estimate [see Fig. [T]. We believe that a physical explanation 
for the scaling of Afe is rooted in the fact that with increasing 
k the PSS \ipk) are increasingly more non-Gaussian, hence 
there is more information not retrievable from second mo- 
ments only. This argument can be made quantitative by eval- 
uating the entropic non-Gaussianity lfl6ll Tfc of \ipk}> which 
simply amounts in this case to the Von Neumann entropy of 
the associated Gaussian state with CM erfe. We obtain, for 
r > 0, 

^-i[^i) + ^*(*±^±I±i)] 

For any r, the non-Gaussianity and the absolute error are very 
close to each other, with Tfe > Afe, and exhibit the same 
scaling with k (see Fig. [2j. This fascinating connection adds 
insight to our analysis and leads us to sum up the results 
achieved so far as follows. Entanglement in ideal photon- 
subtracted states can be measured from the covariance ma- 
trix up to a narrow error that scales with the states' non- 
Gaussianity. 

IV. GENERALIZATION TO MIXED STATES AND 
FURTHER REMARKS 

In most practical implementations, efficient photon- 
number-resolving detectors are not available and the condi- 
tional generation of PSS is achieved by means of "on/off" 
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FIG. 2: (color online). Scaling as a function of k of the asymptotic 
absolute error A™ ax on the estimate of entanglement via second mo- 
ments (dark bars) and of the asymptotic entropic non-Gaussianity 
-j-max Qjgjjj bars) for pure fc-photon-subtracted states \ij)k). 



type detectors which are only able to discriminate the vac- 
uum from a bunch of an undefined number of photons ifTHl . 
This means that a more appropriate description of a class 
of PSS must be in terms of statistical mixtures of the form 
Bk = Y.kPk\i>k){i>k\, where we can define k = J2kPk k as 
the 'average' number of photons subtracted per beam M26I1 . 
For these mixed states not even the logarithmic negativity is 
available in closed form (and its numerical evaluation on a 
computer requires several days for any given squeezing de- 
gree lUlll ). let alone the entanglement of formation. Remark- 
ably, our bounds can be immediately extended to pinpoint 
the entanglement of formation of mixed symmetric PSS from 
the sole knowledge of second moments. We first observe 
that (having zero first moments) the CM transforms linearly: 
o-fe = J2kP k<7k > nen ce it can be computed analytically for 
any probability distribution {pk} from Eqs. ( 1415b . The Ef of 
the corresponding Gaussian state with CM crj., according to 
the extremality theorem fl5tl . stands as a lower bound for the 
Ef of the non-Gaussian mixed PSS qj, ||27jl . On the other 
hand, denoting by 7W the (1, 3) element of the mixed-state 
CM erfe, = J2k Pkl^ k \ the upper bound Eq. (0 is imme- 
diately extended to the mixed case: log[l + 27^ fc '] > EF^Qk). 
The proof follows from the concavity of the log function, the 
convexity of the entanglement of formation, and obviously 
Theorem 1. Namely, log [I + 27W] > Y. k P^°d l + ' 2 l {k) ] > 
J2k PkE{ipk) > Ep(Qk)- The behavior of the bounds in such 
more realistic conditions is shown in Fig.[TJ'(f) for an instance 
with k = 1. We observe, in general, that for reasonable 
modeling of the mixture (e.g., pk following a binomial dis- 
tribution) the loosening of the bounds compared to the ideal 
cases with k = [k\ is negligible: our scheme is efficient 
and robust against the specific source of imperfection con- 
sidered here (the usage of non-phon-number-resolving detec- 
tors). We plan to deepen our investigation in the future follow- 
ing the experimental progresses in the generation of (generally 
non-symmetric) PSS states 12611 . thus properly modeling other 
sources of imperfections (e.g. mismatches or dark counts in 
the photon conditioning) that arise in practical demonstrations 
ll30ll . in order to test the robustness and reliability of our tech- 
niques for estimating entanglement in fully realistic situations. 

In this context, let us briefly comment on the direct imple- 
mentation of our results in experiments. Once a two-mode 
PSS is prepared, one needs to measure the full CM by homo- 
dyne detections as in lfl7l [l9ll . transform it in standard form 
(i.e. extract the symplectic invariants en, 0,2, 7 X , j p ) 1H l20Tl . 
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and then (upon verification that the standard-form CM has the 
structure predicted here: a benchmark for the state engineer- 
ing) readily evaluate our lower and upper bounds to ensure 
an accurate estimate of the entanglement of formation of the 
produced non-Gaussian state. 

V. CONCLUSION 

In this paper, in the spirit of lfl5i l28ll . we have gone be- 
yond the conventional belief that out of Gaussian states the 
covariance matrix plays a marginal role in CV entanglement 
quantification. On the contrary, we demonstrated that clever 
exploitation of such an easily accessible component of the 
state, bears extremely useful and precise information on the 
quantification of non-Gaussian entanglement produced in ex- 



periments, specifically for the important class of (realistic) 
photon-subtracted states ifTol [Till . This is the start of a pro- 
gram which will continue with the systematic investigation of 
quantitative entanglement witnesses for other classes of non- 
Gaussian states in terms of low-order moments. We hope with 
our result to stimulate advances in the engineering and char- 
acterization of non-Gaussian resources, and their exploitation 
for demonstrations impossible to achieve with Gaussian states 
only, in order to explore the actual limits that quantum me- 
chanics poses on the access and manipulation of information. 
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